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Numerical simulations of dynamical gluinos in SU (3) Yang-Mills theory: 
first results. 
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In a numerical Monte Carlo simulation of SU (3) Yang-Mills theory with dynamical gluinos we have investigated 
the behaviour of the expectation value of the scalar and pseudoscalar gluino condensates in order to determine 
the phase structure. Preliminary results are presented as a function of the hopping parameter. 



1. INTRODUCTION 

In the last years there has been a great progress 
in the understanding of the non-perturbative 
properties of supersymmetric gauge theories. Be- 
cause of their highly symmetric nature, super- 
symmetric quantum field theories are best suited 
for analytical studies, which sometimes lead to 
exact solutions The basic assumption about 
the non-perturbative dynamics of supersymmet- 
ric Yang-Mills (SYM) theory is that there is con- 
finement and spontaneous chiral symmetry break- 
ingi- 

1.1. Supersymmetric Yang-Mills theory 

Since local gauge symmetries play a very im- 
portant role in nature, there is a particular in- 
terest in supersymmetric gauge theories. The 
simplest examples are SYM theories, which are 
supersymmetric extensions of pure gauge theo- 
ries. We shall pay our attention to the SYM 
action with iV = 1, where N is the number of 
pairs of supersymmetry generators Qid , Qia (* — 
1, 2, . . . , TV). This theory is a Yang-Mills theory 
with a Majorana fermion in the adjoint represen- 
tation. 

The action for such a = 1 SYM theory with 



a SU {Nc) gauge group is given by 

^ = F%{^) + \'^\x)^,V,^%x), (1) 

where '^"■{x) is the spinor field, and F^^(a;) the 
field strength tensor, a G {1, . . . , — 1}. 

Introducing a non-zero gluino mass rrig breaks 
supersymmetry "softly" . Such a mass term is 



(2) 



Here in the first form the Majorana- Weyl compo- 
nents A, A are used, in the second form the Dirac- 
Majorana field 4". The Yang-Mills theory of a 
Majorana fermion in the adjoint representation 
is similar to QCD: besides the special Majorana- 
feature the only difference is that the fermion is 
in the adjoint representation and not in the fun- 
damental one. As there is only a single Majorana 
adjoint "flavour" , the global chiral symmetry of 
A^ = 1 SYM is U{l)x. 

The t/(l)A-symmetry is anomalous: for the cor- 
responding axial current ~ ^'7^75^', with a 
gauge group SU{Nc), we have 
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However the anomaly leaves a Z2Nc unbroken: 
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this can be seen by noting that the transforma- 
tions 



are equivalent to 

nig m^e^^"''^^ 
and 
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(4) 
(5) 
(6) 



where 6sym is the ^-parameter of the gauge dy- 
namics. Since ©sym is periodic with period 2tt, 
for rrig = the U{l)x symmetry is unbroken if 

^ = <^fc = ^, (fc = 0,l,---,2A^c-l) . (7) 

The discrete global chiral symmetry Z2Na is ex- 
pected to be spontaneously broken to Z2 by the 
non-zero gluino condensate (x)\E'(x)) 7^0. The 
consequence of this spontaneous chiral symmetry 
breaking is the existence of a first order phase 
transition at zero gluino mass rrig = 0. In the 
case of A^c = 2, there exist two degenerate ground 
states with opposite signs of the gluino conden- 
sate. An interesting point is the dependence of 
the phase structure on the gauge group: instan- 
ton calculations |^ at 8sym = give Nc degen- 
erate vacua {k = 0, ■ ■ ■ , Nc — I) with 
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The coexistence of Nc vacua implies a first or- 
der phase transition at nig = 0. Recently Kovner 
and Shifman have suggested the existence of an 
additional massless phase with no chiral symme- 
try breaking |Q. In the case of SU{3), there are 
at least three degenerate vacua and for nig < 
we expect that 8sym = tt. 

2. LATTICE FORMULATION 

No lattice gauge theory exists with an exact su- 
persymmetry. This is because lacking lattice gen- 
erators of the Poincare group, it is impossible to 
fulfill the (continuum) algebra of SUSY transfor- 
mations. Another problem is represented by the 
balancing between bosonic and fermionic modes 
required by SUSY: the naive lattice fermion for- 
mulation produces too many fermions. 



Curci and Veneziano ||5|] have proposed a sim- 
ple solution: instead of trying to have an exact 
version of SUSY on the lattice, the requirement 
is that, like chiral symmetry, it should only be 
recovered in the continuum limit, tuning the bare 
parameters (gauge coupling g, gluino mass nig) 
to the supersymmetric point. 

2.1. Actions 

The Curci- Veneziano action of = 1 SYM is 
based on Wilson fermions. The effective action 
obtained after integrating the gluino field is given 
by 
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Scv = /3^(1 - -TrC/pO ~ - log detQ[U] . (9) 



The fermion matrix for the gluino Q is 



K ^ (5y^a; + /i(l+7^)K„, 2,^(10) 



with the gauge link in the adjoint representation 
Vvu,xi^i[U] = Tr(C/t^Ti,C/j;^T„). 

2.2. Monte Carlo simulation 

The renormalized gluino mass is obtained from 
the hopping parameter K as 
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ZmiafJ,)nio~g . (11) 



Here Kq = Kq{(3) gives the /3— dependent posi- 
tion of the phase transition and /x is the renor- 
malization scale. The renormalized gluino con- 
densate is obtained by additive and multiplicative 
renormalizations: 

Z{a^i) [< l'(x)'I'(x) > -6o(aAi)] . 

A first order phase transition should show up 
as a jump in the expectation value of the gluino 
condensate at i^T = Kq. By tuning the hopping 
parameter K to Kq for a fixed gauge coupling (3 
one expects to see a two peak structure in the dis- 
tribution of the gluino condensate. By increasing 
the volume the tunneling between the two ground 
states becomes less and less probable and at some 
point practically impossible. It is possible to see 
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this phase diagram in our simulations by measur- 
ing the chiral and pseudo chiral gluino conden- 
sate: the order parameter of the supersymmetry 
phase transition at zero gluino mass is the value 
of the gluino condensate 



(12) 



Additionally, for K > KQ{mg < 0) a spontaneous 
CP-violation, indicated by a nonvanishing pseudo 
condensate < "if {x)j5'i/ (x) >^ 0, is expected. 

We determine the value of p on a gauge config- 
uration by stochastic estimators 



(13) 
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Outside the phase transition region the observed 
distribution of p can be fitted well by a single 
Gaussian but in the transition region a good fit 
can only be obtained with two Gaussians. For 
SU{2) results are shown in 

The hopping parameter Kq, corresponding to 
zero gluino mass, is indicated by a first order 
phase transition which is due to the spontaneous 
discrete chiral symmetry breaking Zg — > Z2. 



Figure 1. The preliminary results for the gluino 
condensate. 

We have investigated the dependence of the dis- 
tribution of the gluino condensate and the pseudo 
condensate as a function of the hopping parame- 
ter, starting from a lattice volume ■ T — ■ 8. 
This lattice is, however, still not very large in 
physical units. Therefore the expected two-peak 
structure is not yet very well developed, never- 
theless we have high statistics. For K = 0.195 



fig. H shows the distribution of the gluino con- 
densate. The distribution indicates that we are 
near the phase transition. Outside this region, 
we can fit the distribution with a single Gaus- 
sian. Presently, we are calculating on a bigger 
lattice volume (L^ • T = 6^ • 12) in order to sepa- 
rate the two-peak structure. On the other hand, 
our present results on the smaller lattice do not 
show any signal for a pseudo condensate. 
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Figure 2. The probability distribution of the 
gluino condensate for K = 0.195 at /3 — 5.6 on a 
43 • 8 lattice. 
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